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Abstract 

We study the statistics of a system of N random levels with integer val- 
ues, in the presence of a logarithmic repulsive potential of Dyson type. This 
probleme arises in sums over representations (Young tableaux) of GL(N) in 
various matrix problems and in the study of statistics of partitions for the 
permutation group. The model is generalized to include an external source 
and its correlators are found in closed form for any N. We reproduce the 
density of levels in the large N and double scaling limits and the universal 
correlation functions in Dyson's short-distance scaling limit. We also study 
the statistics of small levels. 
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1 Introduction and definition of the generalized 
model 



The theory of random matrices leads often to consider character expansions, i.e. sums 
over irreducible representations of groups such as GL(N) or U(N) |1|, |2], [| |, ||, |6], [7|, 
§, [|. Similar sums play also a role in probability theory in particular in recent 



studies of the distribution of cycles in the group of permutations [jlTJ, [12|, [13], [14], [L5 
This has led to the study of probability distribution functions (PDF) for N non- 
negative integers hi, ■ ■ ■ , h N of the form : 

V a (h lr --,h N )~A 2 (h)a£>ti h >> (1) 

where A(h) = Y\ m> j{h m — hj), and < a < 1 is a real parameter. 

Let us generalize this distribution to an ensemble characterized by N parameters 
«i, • • • , «tv, in the following way : 

'P[a 1 ,-,a N ](hi, • • • , h N ) = ^A(h) Xh ({a}) (2) 

where XhM = ^ ^4 °. k is the Weyl character of a diagonal GL(N) group element 
A(a) 

ax, • • • , a^v of a given irreducible representation i£. i? is fixed in terms of the highest 
weights rrik = hk + k — N, k — 1, • • • , N. (In principle the corresponding h^s are 
strictly ordered, but in view of the symmetry of the weight (|1|) this restriction may 
be ignored in the sums). The constant Z is defined by the normalization condition 

oo oo 

A(/0xfc({a}) (3) 

/il=0 hjy=0 

Note that in the limit of coinciding a's — > a the (non-normalized) distribution 
(0) reduces to (JTJ) ; this follows from the well known formula: 

tti(a) - ft « E ^ mfc dim {m} = *W a E>^-™0 (4) 

fc=l 1 lfc=0 fc - 

where dini{ m } is the dimension of a representation given by the highest weights 
m 1 ,---,m N . 



There are several reasons for this generalization. First it leads to simple exact 
formulae, as in the case of random matrices coupled to an external matrix source 



16| , [I7| ; the N parameters c^'s play here the role of the eigenvalues of the source. 
Here also, even when the source is a simple multiple of the identity, meaning now 
that all the a^s are equal to a single a, the final formulae are explicit and simple. 
Furthermore these parameters a^'s provide a powerful check of universality. Indeed 
it is found here that at generic points, at which the density of levels is non-singular, 
the correlations are, in the proper Dyson scaling limit, insensitive to the specific 
probability distribution. Singular points fall also into universal classes, as for the 



usual matrix models Varying the external parameters one can indeed tune 



various singular classes : examples are given in the subsequent sections. The machin- 



ery developed for integrable systems, and used in |TT|, [L2|, [13], [14], [L5[] for the single- 
parameter model (|l|), ought to be be applicable to our generalized model (0) as well. 

Finally the distribution probability ([!]) provides a natural measure on the Soo 
group of arbitrary permutations ; the m^'s with k = 1,2, •• -,N, are the lengths 
of the cycles of a permutation class consisting exactly of N cycles. One can also 
interpret this distribution probability as defined on the (infinite) set of all Young 
tableaux with rrik boxes in the fc'th row. 

In our case, (0) is a natural multi-parametric generalization of ([I]), which may 
now be interpreted as a specific coloring of Young tableaux. The boxes of the Young 
tableau characterizing a permutation, are colored in N colors in a way which is 
explained below ; the k-th color is weighted with a^. For instance, for N = 2, 
there are two rows of lengths (highest weights) m x > in the corresponding Young 
tableau characterizing a class of permutations with 2 cycles, and for the character we 
have the following finite sum of positive terms: 

/ \ x-^m±— mo mi— k mo+k 

= a^af 2 + + • • • + «r +1 «r _1 + «r«r (5) 

This may be interpreted as a coloring of a tableau, in which the first ni\ — k boxes 



of the upper row have a color 1 and all the other boxes have color 2 ; we sum over k 
with a factor af 1 af 2 where #1, #2 are the numbers of boxes ("areas") of the Young 
tableau of a given color. 

To generalize this formula to any N, we have to expand a general character into 
a sum of monomials in the a's (generators of the maximal torus of GL(N)): 

N 

X {m} = E "(Oil («*)'* ( 6 ) 

W{m} k=1 

where the finite sum over the positive integers {h, ■ ■ ■ ,In) characterizing the ele- 
ments of the Lie algebra of the maximal torus, or the weights of representation R, is 
restricted in a specific way by the shape of the corresponding Young tableau of the 
representation R(mi, ■ ■ ■ , m N ). The positive integers n{l} are called the multiplicities 
of those weights. 

To make all this more explicit, and to give it a nice probabilistic interpretation, 
let us introduce the Gelfand-Tseytlin scheme (GT scheme), which provides an or- 
thonormal basis of states for a given representation R of GL(N) characterized by 
the highest weights m,, % = 1, • • • , N. Every state vector of R is characterized by 
N(N — l)/2 positive integers m^j, k = 1, ■ ■ ■ , N, j = 1, ■ ■ ■ , k. The first ones, the 
m^r j = rrii, i = 1, • • •, N, are simply the highest weights. The subsequent ones are 
given integers, restricted by the inequalities: 

m^i > rrij-ij > m jti+1 , i = 1, • • • ,j (7) 

The basis vectors of this representation are thus characterized by a triangular array 
usually depicted as |TB[ 

m N}1 m N)2 ■ ■ ■ m N)N 

Wjv-i,i ■ • ■ mjv-i.iv-i 

m N<N (8) 
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in which any m in a given row lies in-between the two m's next to it in the previous 
row. We will use the definition of a character as a trace of a group element in its 
diagonal form, for a given representation R: Xn({ a }) — TrRFIjli^j JJ > where Tjj is 
a diagonal generator in the Lie algebra of this representation. In the GT basis the 
values of T$ are given as (see for example |18| ) 

3 3-1 

T jj = l i = J2 m 3i ~ £ m i-l,i ( 9 ) 
i=l i=l 

The expansion (Q) is now explicitely given by the formula: 

AT 

x*(M)= e n«y {m} ( io ) 

{GT} B j=1 

where the sum is taken over all GT schemes (states), satisfying the restrictions ([?]). 
Of course in this sum the same monomial may appear several times. The number of 
times each monomial enters in the sum is the multiplicity n{l} of eq. (|6]). 

This formula may be given an interesting probabilistic geometrical interpretation, 
in terms of sums over the colorings of Young tableaux. Let us define a coloring in the 
following way (see figure 1 for an illustration of the N = 3 case). The color number 
one is given to the last mjv,i — w&jv-i,i boxes of the first row, the last mjy,i ~ wijv-i,i 
boxes of the second row, etc., the m^^N boxes of the last row. Color number 2 is 
given to the boxes numbered mjv-i.i + 1 to m N _ 2 ,i in the first row, m N _i i2 + 1 to 
mN -2,2 in the second row, etc., the first 77ijv-i,jv-i boxes of the row before last ; etc., 
the iV-th color is given to the first m^i boxes of the first row. At the end, the Young 
tableau is a "brick wall" made out of colored bricks, each brick at an upper level 
relying on zero, one or two lower bricks of lower color indices, as depicted on fig.l. 
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m 33 = 2 




Fig. 1 : Coloring of the Young tableau for the term in the distribution probability proportional 

h h h 
to a 1 a 2 0C3 J . 



Note that if all the a&'s reduce to = 1, one is simply summing in (10) over 
each state with weight one ; therefore one recovers the dimension of the representation 
dirriR as stated in ([J). 

If one substitutes the expression (0) for the character into @, one obtains a 
probability for the colored Young tableaux (related to permutations of colored ob- 
jects). 

In the next sections, we derive explicit (for any N) formulae for the correlators of 
the distribution (@) and study various scaling limits. 



2 Density of discrete levels 



The density of levels is defined as 

in which the brackets denote the averaging with respect to the probability distribu- 
tion (Q). The variable Na is a priori an integer but, in the large iV-limit, a will be 
considered as a continuous variable. Let us calculate Ui(t) = (-^J2k^ hk )h by per- 
forming explicitly the sums over the h's. Using the antisymmetry of A(h) in fl2|) we 
can rewrite it as 

hi 



£M*) = c££<**%^r = ^EEHTII («2>)*" "™ < 12 > 

k h ^V*) ^V*) k P h rn 



where the determinant is represented as a sum over permutations P of the levels 
hi, ■ • ■ ,h n and by definition 

«W = t 5 ^a m . (13) 

The sums over h's are now independent and can be calculated using the formula : 

00 1 1 1 / 1 \ 

E h p a h = (ad a y- = - Q p 

^ 1 — a p! 1 — a VI — a) 

where Q p (x) is a polynomial of degree p, whose coefficient of highest degree is one. 

Since for any set of such polynomials det Pi9 Q p -i(x q ) = A(xi, • • • , xn) we obtain from 



(0): 



C ^ A ((1 - qW)- 1 ) i 



A ^ r n m f i - N k 



n AT / 



1 — tai 



a k - aij 



where prime in the product means that the term k = j is omitted. The constant C 
has been determined, through the normalization condition Z7i(l) = 1. 

This quantity has an elegant contour integral representation, similar to the one 
found in [16j for integrals over hermitian random matrices with an external matrix 



source. Indeed, the integral 



TT . . 1 r du ( 1 — u\ N \-rtu — a-j 



iV(t — 1) / 2z7ru \1 — tu) u — a 
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over a contour in the complex u-plane which surrounds the N poles a/s (and not the 
origine), reproduces the sum fll4l) . For the density of levels (0), since 

r dt Uitf) 

P(a) = /w^' (16) 

where the integration contour encircles the vicinity of the origin, we obtain, using 
If dt r du ( 1 — u \ N -f-f tu — a* 

p{a) = Nf^t^{t-i)f^u{—u) n^-- ( i7 ) 

1 N 1 

Similarly for the resolvent, defined as G(z) — (— — — -), we can use 

y J N ~r[z — h k /N 



POO 

G(z) = / dre~ TZ tfi(e T/JV ) 
Jo 



and obtain the integral representation: 

rM _ 1 f°° dr _ T r du f l-u \ N u e T / N u - Oj 
[Z) ~NJo (e-/ N -lf J 2i7ru\l-e^uJ A / u - a j 1 j 

After the change t — > t/u in fllTD we get (er = p/N): 

; iV 7 2i7rtP +1 - 1 / (t - 1) 7 2i7r Y (u-aj)t-u V ; 

or, inflating the contour of integration in u and changing u —>■ 1/w,: 

, „ n If dt t-t (a* — t) f dw -n- (1 — w) 1 

Expanding the last factor we get a finite sum representation for the correlator: 

p(p) = ^t^N k (^ P N k+ \<y) (22) 

iV k=0 

where the polynomials in a's C 9 N (a) and £/v-( a ) are defined as 
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The contours for both ( p3| , f2~4| ) encircle the pole at x = 0. 



Note that this expression is exact, for any N and any set of a's. In particular 
for the density, in the simplest case a,\ = a 2 = ■ ■ ■ = ot^ = &, most studied in the 
literature, this leads to a simple integral: 

. /Ar . -1 r dt (t - a) N r duvP +1 (u - l) N 1 

or, after inflating the contour and changing u — > a/w 

, .. T . r du r dt f (1 — u)(a — tu)\ 

= -f^-ut 2i^(t -1) ( (1- tu)(u - a) ) ■ (26) 

The polynomial (p3|) reduces now to: 

^)=f^-r(^;) N (27) 



2irix q+1 \l-x 
and the two sets of polynomials are related here simply by 

C%(a) = a N -«C q N {a). (28) 

In the next sections we will study the large N limit of this density and, for the 
case of equal a's, compare it with a direct computation based on the solution of a 
Riemann-Hilbert problem given in the appendix. 

3 Pair correlator of discrete levels 

The simplest interesting correlator, which is conjectured to obey a short distance 
universality for nearby levels is the pair correlator: 

U 2 (t h t 2 ) = j^(J2t h i k t2 l ) (29) 

k^l 

or rather its connected part: 

K 2 (t u t 2 ) = U 2 (t u t 2 ) - U x (t x )Ui (h) (30) 
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A similar calculation yields 

TT (. , v 1 v-> (hoik - t 2 aii)(a k - on) 

iV 2 ^ (ti« fe - - t 2 ai) 

-r-r (ftm ~ glOfc) rr (ttm ~ ^2«z) 

m ^ fc (a ro -a fc ) (am-aj) 

\ L \l-a k tJ \l-a k tj V ; 

from which one derives again a contour integral representation over two complex 
variables (similar to the representation for the hermitian matrix model found in |T7|]): 

1 r r dudv ( 1 — u \ / 1 — v x A 



(2iix) 2 (u — t 2 v)(v — t\u) \\ — t\uJ \1 — t 2 v 



tt {hu - aj){t 2 v - aj) 
Finally, for the density correlator 

pip, q) = (E <W* E <W> = i i . d li„.\ ^' t2 ) ( 33 ) 

we obtain ( after changing variables from ti to t\ju and t 2 to t 2 /v), a factorized 
formula, again in similarity with the hermitian matrix model, ||17||: 



pip, q) = -R(p, q)R(q, P) (34) 

where 

, > f f dudtt- p u q - p - 1 /u-l\ N ' (t- o,) , s 

«M = ff (2,).(i- t ) (— ) (35) 

or, going back to the previous non-scaled variables: 

. , f f dudtt~ p u q T-rit — aA (u — 1) 

*«»=-//wMnFli4 (36) 

Another useful form of R(p, q) can be obtained (similarly to (pTj) ) by the contour 
of integration in u in ( |36l) and changing w — > 1/w to catch the poles at w = 0: 

/■ f dwdtt~ p u~ g (aj - 1) (1 — to) 1 

fl <^> = // ^>^n t ; ( ; _ ^ (1 _ tw) (37) 
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Expanding the last factor we get a finite sum representation for the correlator: 

inf(p,q) 

R(p,q) = ar N+ * £ C p N ~ k (a)C%-\a) (3f 

fc=0 



where ^(a) and £^r(a) are given by and 

In the next sections we will study these formulae in the large N limit, for the case 
of large Young tableaux. 

4 Study of the density in the large N limit 



In the large N limit the formula for the resolvent fll9|) gives 

G M = f t/^H'-iV ^!) (39) 

where G (u) = J2j=x is the resolved of distribution of parameters acj of our 
problem. Integration over r leads to the eq.: 

G ^=ftuH l "-T^- uG ^ )] ) (40) 

Differentiating the last eq. in a, taking into account the contribution of the pole 



at 



ii 

a = + uG (u) (41) 

1 — u 



and integrating back in a we obtain an explicit equation for the resolvent 

G{s) = \nu(s) (42) 



where u(s) is a solution of the eq. ( 41 ) . The constant of integration in ( 42]) is chosen 



to be zero in order to match the assymptotics a — » oo (u — > 1) 



If we eliminate u from the eqs. (|4lD and fl42|) we obtain the following functional 
equation for G(a) 

e G{a) = 1 + - + -G (e G(<7) ) (43) 
o a 
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This functional formula is very similar to those found by the direct analysis of the 
saddle point equations in || for the heat kernel on the group SU(N) in the large 
N limit (see also || [7] for many other similar results). 

For the particular case a\ = a 2 = ■ ■ ■ = «jv = a we have G (u) = and G(a) 
is given by 



G{a) = log I -!- I (a + l)a - (1 - a) - (1 - a) 



1 - V« ir 1 + y/a, 
\a - - =\\a 



' (44) 
Since on the real axis 

G(a ± 0) = ReG(a) ± iitp(a) (45) 

we see from ([PJ) that 

p(a) = 1, for < o < b (46) 



p(cr) = log y(a + l)cr — (1 — a) — (1 — «)y (cr — 6) (a — a)j j ,for b < a < a 

(47) 

where 

a = =, o = —. 48 

This coincides with the direct saddle point calculation of the Appendix and repro- 
duces an old result of A. Vershik and S. Kerov [0 on the limiting shape of the Young 
tableau. 

Another interesting problem corresponds to the statistics of the weights close to 
the upper or lower end of the distribution |L^, 0] etc. studied in the next section. 

5 Double scaling limits 

The density function contains three special points ("end points"): z c = 0,a,b, near 
which it exhibits a universal behavior. We will study the close vicinity of the points 
a, b such that the deviation Az from one of these points scales with iV according 
to the power law Az ~ jV~ 2 / 3 typical of the double scaling limit for the generic 
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distributions of the eigenvalues of matrix models [B5[. The scaling for a — > 1 in the 
vicinity z ~ 6 ~ will be different and will be considered afterwards. 

It is useful to put in ( |I9|) t = e T/W = 1 + 9/N. This leads to the integral repre- 
sentation: 

nt \ 1 f°° de I du f-, u °\ N -nf, U °\ /a^ 



Exponentiating the integrand and expanding in 1/iV we get: 

G(z) = f eoTW) / 5; exp 1 ~ " ~ uG " (u] 

+ w (-■ + jt^w + u2G '° (u) ) - + ^{- z+ ir^ + r' G ^) + ■ ■ ■ i < 50 > 

The end points a, b are defined by the equation ([0]). The singular behavior arises 
when the linear in 9 term in ( |50]) develops a double zero at the the end point (i.e., 
when its derivative in u is zero): 

G (u) + uG' (u) + 1 = (51) 
(1 — u) A 



One sees immediately that if the conditions ([41) and ([51]) are satisfied the quadratic 
in 9 term is also zero. In the vicinity of the end points z = a or z = b we can expand 
in the exponential in powers of A = z — a (or respectively (z — b)) and Su = u — u c : 

- L 00 iB f exp i - ( A -\ A ^ 2 ) + £ < A + ^ + + ■ ■ ■ i < 52 > 

with A = £ + - ^— r + M c Gn. In the double scaling regime, when A ~ iV~ 2 / 3 , 

uj (1 - u c ) A 

the quadratic term in 9 2 A/N terms is negligible, and after Gaussian integration in u 
along the contour of the stationary phase, we obtain an Airy-like function of AiV 2 / 3 

& nt \ f°° d9 -AN 2 / S e-B8 3 



This gives for the density a double scaling expression: 

p(cr) ~ A 1 / 2 /(AAT 2 / 3 ) (54) 
12 



where the function f(x) of the double scaling parameter can be defined nonpertur- 
batively by an appropriate change of the integration contour in (|53|). 

For the last end point of the distribution, namely the vicinity of z = 0, in the 
simple large N limit p(z) = 1 around this point. But it will be not so any more in 
the special double scaling limit in which a approaches 1 and the interval (0,6), on 
which p remains fixed to one, shrinks to zero. We will consider this special limit in 
section 7. 

6 Study of the density correlator for nearby levels 

Now we shall study the quantity R(p, q) in the large N limit for a separation of p, q of 
order one : \p — q\ — 1, 2, 3, • • •. The result is expected to be universal (if we express 
p — q in terms of the local level spacing p{p)~ 1 )- 

The study of the large N limit of R(p, q) is very similar to the study of the density 
p(p) in the previous section. Namely we put again t = e T l N and retain in fl36|) only 
the terms of the order ~ 1. The integration over r yields the following large N limit 

r(p, t ) = Kfe/iv)]'"' - M./AQ]"-" (55) 

\p-q\ 

where u±(a) = e G ^ are two conjugated solutions of (|4lD corresponding to two choices 
of the sign in Q35| ) when approaching the real cx-axis. Hence we get on the real axis: 

R(p, q) = 2ie^W ^Mg/^))P-gl) (56) 

\p-q\ 

The level correlation function is thus 



p[M) ^MilN)b-i\) (57) 

\p — q\ 

It reproduces the well-known universal formula for the level correlation function with 
Dyson repulsion law (for the correlations of eigenvalues of the unitary ensemble of 
matrices) . It does not change even on the saturated part of the corresponding Young 
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tableau where p(q/N) = 1. The only difference with respect to the hermitian ensem- 
ble is that in our case of discrete levels is that the this correlation function makes 
sense only at discrete values of the distance between levels \p — q\ = 1, 2, 3, • • •. 



The authors of [|T^, [13], [TJ], [15] came to the same conclusion in a particular case of 
equal a's. Our generalization to any collection of a k s shows once again a remarkable 
universality of the classical result (|57|) . 



7 A special large N limit: small weights in very 
large Young tableaux 

In this section we study a new singular scaling regime, in which all the a k — > 1 and 
N — > oo, so that the parameter p = N(l — — ^ remains finite. 



In this limit the polynomial (a) defined by fl23f) becomes 



N k 



£n(<*) -> * e P/V-V = e -P <k ^ ' e~ pz = e- p MJp), (58) 

N K J J 2rritP+ 1 J 2mzP +1 pKHJ K J 

with 

M p {p) = L p (p) - L^p) (59) 

in which L n (p) is the standard Laguerre polynomial of order n, normalized to L n (0) = 
1. Similarly, in this limit, 

C%\a)^M p {p) (60) 

The formula fl22|) becomes: 

Vooip) = e"' £ M k (p)M k+1 (p). (61) 

fc=0 

In the limit p — > of ultra large Young tableaux, we have M p (p) ~ — p and the 
formula flBTD gives explicitly 

^oo(p)=p-(p+l)p 2 + 0(p 3 ) (62) 
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In full analogy with these formulae for the density we can deduce from 
correlation function of the small weights p and q for large Young tableaux: 

^oo(p,g)= E M p _ k {p)M q _ k ( P ) 

fc=0 

which gives for the limit p — > of ultra large Young tableaux 

Kodp, q) = -p + inf(p, q)p 2 + 0(p 3 ) 

8 Conclusion 

In this work we have studied discrete statistics of the kind Such character ex- 
pansions (CE) appear in many situations in which one sums over the irreducible 
representations of GL(N) : the hermitian matrix models, circular ensembles, in vari- 
ous partition functions and loop averages of gauge theories , such as those encountered 
in two-dimensional quantum chromodynamics |]T], 0. The systematic large N saddle 
point analysis of various CE's was first proposed in || in relation to the calculation of 
partition functions and Wilson loops of QCD2 on the sphere, and developed further 
in j|, |5|, 0. It appeared to be the only effective method in the study of the combi- 
natorics of the so called dually weighted planar graphs (DWG) |J first introduced in 

I- 

Although we studied only a particular example of such models our results suggest 
that the universal properties of various correlation functions observed for our model 
will hold for many other models defined as multiple sums over discrete levels (random 
ensembles on the Young tableaux). 

The explicit formulae, exact for finite N, which have been obtained here-above 
for the density of levels and for the correlation functions, provide a powerful handle 
on the study of various universal scaling limits. 

It would be nice to apply our methods to more complicated characteristics of the 
model, such as the probability of distribution of the length of various rows of the 
Young tableau, (which in the case of equal a's shows, according to the authors of 
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H) the 
(63) 

(64) 



IT I, a remarkable relation to the Painleve II equation). Unfortunately, for such a 
quantity, we have not succeeded to write it in the form of a finite contour integral 
representation, similar to that used in our paper for the correlation functions. 
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9 Appendix : Density of levels by a direct saddle 
point calculation 

It is interesting to compare the result obtained in section 5 with a direct calculation in 
the large iV-limit. For the simplest model in which all the a's are equal, the method is 
straightforward ; it is based, as usual, on the solution of a Riemann-Hilbert problem. 

The probability weight of a given set of random Tt's, is proportional by definition 
to a^ hk A(h u ---,h N ) . Is we write it in terms of the density distribution 

1 N 

p(x) = -^5(x-h k /N), (65) 
this weight is proportional to 

-p * + /**^ - m) ) . x- « * *. ^ 

tary ensembles of random matrices, with their characteristic logarithmic repulsion, 
is manifest. In the large- N limit the weight is maximum for a density which satisfies 
the condition 

log a + IP f dy-^- = (66) 
J x — y 

. This condition holds on the support of the measure, which is yet to be determined. 
We shall make the ansatz (proposed in 0] in a similar situation) that p(x) remains 
equal to one for < x < b, is some function of x in the interval b < x < a, vanishes 
at x = a and remains zero for x > a. The consistency of these conditions will be 
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checked later. The equation for p{x) on the interval b, a becomes 

log + log + P r dy^- = (67) 

\x — 0/ Jb x — y 



We introduce the resolvent 

ra ^ pjy) 

z-y 

, and 



G M = r^ < 68 > 



fl w=r*^ (69> 

. The function H(z) is analytic in the cut-plane from b to a. On the cut it satisfies 

ReH(x) + log ^| = 0. (70) 
One verifies readily that the function 

is the only fuction which satisfies the analycity requirements and the behaviour at 
infinity. If we consider a large circle in the complex u-plane, the behavior at infinity 
implies that 



' togPqUo. (72) 



z- u y/(u-b)(u-a) \ u - b 
Shrinking now the contour, and collecting the various singularities, we obtain 

The integration is then elementary and we end up with 

1 Jb(z — a) + Ja(z — b) 

G(z) = - -log a - 2 log V-^ jLJU L. (74) 

Z \Jz{a-b) 

The behavior at infinity G(z) ~ - fixes a and b ; one finds 

a = i±^and b = (75) 
1 - V" 1 + v 7 " 
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One finds then 



G(z) = log— ^ f(l + a)z - (1 - a) - (1 - a)Jz - a)(z - b)\ 
2az \ v / 



(76) 



in agreement with the result of section 5. The density of levels is then obtained as 



p[x) = ImG(x + iO), 



7T 



(77) 



and one verifies our assumptions : p(x) = 1 for < x < b, vanishes for x > a and in 
the interval (6, a) 

t \ 2 b(a — x) . s 

pix) = — arctan , — — . (78) 

7r \ a(x — b) 

This result has been first obtained by A. Vershik 
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